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Measurements of shear stress for various shearing How histories at 25 °C are correlated through 
the BKZ elastic fluid theory. The data are on 10 percent P1B solution in cetane. The histories include 
single step stress relaxation, suddenly applied steady shear, steady shearing flow, stress relaxation 
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1. Introduction 

In this paper we present the results of a series of 
experiments on a single sample of polyisobutylene 
solution subjected to various histories of simple shear. 
The work was done as a part of a larger program to 
study the appropriateness of a fairly general constitu- 
tive equation, the Bernstein, Kearsley, Zapas (BKZ) 
elastic fluid [l], 1 in describing a wide range of non- 
linear viscoelastic behavior for a range of materials. 
This larger purpose influenced both the choice of test 
material and the choice of experiments. In the first 
place, the studies were confined to a single sample so 
as to provide a test of the correlations between experi- 
ments predicted by the assumed constitutive equation. 
The sample was chosen to show essentially complete 
relaxation within the time scale of the experiment, 
since earlier studies [2] had been made on materials 
with longer relaxation times. The deformations are 
various histories of simple shear, because they are 
experimentally convenient deformations and because 
they are convenient to analyze in terms of the assumed 
constitutive equation. 

Some measurements of normal stresses were made 
on this sample, but we encountered serious discrep- 
ancies when we attempted to correlate them with the 
shear data. This trouble has been traced to the recently 
discovered "pressure hole error" [3] and the difficulty 
resolved, but we have postponed publication of the 
normal stress data to avoid further delay in the 
presentation of the shear stress measurements. 

In the course of these studies some measurement 
errors, customarily ignored, were found to be sur- 
prisingly large and some correction methods were 
worked out and are presented here. A specialized form 
of the BKZ elastic fluid was used in this paper to 



'Figures in brackets indicate the literature references at the end of this paper. 



correlate our experiments with very good success. 
Although this special form is neither unique nor exact, 
it has been chosen to be consistent with a wide range 
of phenomena beyond that described in this paper, 
and should be of considerable practical use. 

2. General Considerations 

The class of motions which we shall consider will be 
limited to simple shearing histories. 

In the incompressible isothermal form of the BKZ 
theory [1, 4| the value of the shearing stress at time t 
is given by the following relation 

<r(t)= (2.1) 



■/: 



W 7 *(y(/)-y(r), t-T)(y(t)-y(T))dr 



where y(r) is the amount of shear at time r and W* 
is a function of strain and time. 2 If W * is known, the 
shearing stress a(t) for any simple shearing flow can 
be calculated through eq (2.1). Of course, the form of 
W* is not specified by the theory and in general the 
function W* must be determined through experiments. 
A good idea of the function W* can be obtained from 
examining a single step stress relaxation history in 
simple shear. In such a deformation the shearing 
history is given by: 



y(r)=o 



for r < 



(2.2) 



y(r) = y= constant for r > 0. 

By substituting £ = t — r and (2.2) in (2. 1 ), we get 

^T = - f " ^*(r, f )4f = W{y, t) (2.3) 



- Equation (2.1) is equivalent to eq (3.11) of Bernstein [4] except tor a change of notation. 
// * of this equation is equal to — 2(£// t + Ut t ) of Bernstein, and our o-(/) is Bernstein's cr n . 
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where 



and 



r*(y,o 



_dW(y, t) 
dt 



(2.4) 



W(y, oo) =0. 

In this paper we shall refer to W as the relaxation 
function. From (2.3) it follows that W or W* may be 
obtained from single step stress relaxation measure- 
ments in shear at different strains and from this the 
stress for any simple shearing flow history may be 
predicted. This approach is similar to the one used 
for simple extension histories [2], 

Bernstein [5] has shown that the function 



%f(y,t) = W(y, *) 1 + 



d In W(y % t) 



d In y 



(2.5) 



may also be obtained from experiments of suddenly 
applied steady shear. Thus, in principle, we should 
be able to obtain W(y, t) from a series of experi- 
ments, and then proceed to check the BKZ theory by 
comparing measured and predicted stress for simple 
shearing histories. 

Of course, one can never realize experimentally 
these shearing histories exactly as assumed. From 
some of the experiments which approximated suddenly 
applied constant rate of shear history, discussed in 
section 6, we were able to get a rough approximation 
of the relaxation function. Applying corrections which 
will be discussed later using an iterative scheme, 
we recalculated a function W which is consistent with 
all our experiments which include measurements of 
viscosity as a function of the rate of shear, stress 
relaxation after shear (for different rates of shear), 
measurement of stress as a function of time for sud- 
denly applied steady shear, and single step stress 
relaxation experiments. These results, we felt, justi- 
fied the use of an expression for W(y, t) which can 
describe all our experiments and which is a special 
form of a more detailed expression consistent with 
the behavior of other materials and other deformations, 
viz: 



W(y,t) = G(t) 



1 



-+- 



1+ 3 y2 



0.145 | y I 

(l+0.06y 2 ) 5 ' 2 



(2.6) 



where G(t) is the shear stress relaxation modulus. 3 
It is evident that there are many other forms which 
are consistent with our experimental data. Our motiva- 
tion for choosing form (2.6) is its relative simplicity 
and its consistency with a strain potential function 
[1J. Furthermore, form (2.6) is useful in process control 
engineering since an immediate translation from the 
linear to the nonlinear behavior can be made as a 



3 Notice that G'(/). the shear relaxation modulus, is the limit of relaxation function 
W(y, t) as y goes to zero. 



function of temperature, molecular weight, and 
concentration. 



3. Experimental 

The data reported in this paper were obtained on a 
10-percent solution of polyisobutylene (vistanex L-100, 
Enjay Chemical Co.) 4 in cetane. A Weissenberg 
Rheogoniometer was used to shear the sample be- 
tween a flat plate (7.5 cm diameter) and a cone such that 
the angle of the gap was 0.0268 radians. The cone was 
at the bottom and was connected to the driving shaft. 
The plate was at the top and was connected to a torsion 
bar which was used to measure the torque. For most 
of the experiments a torsion bar of Vs-in diameter was 
used. The stress-time measurements and the dynamic 
response were recorded on an oscillograph. 

The chamber enclosing the cone-plate assembly 
was kept at a temperature of 25.0 °C ± 0. 1 °C. Ambient 
temperature was controlled at 25.0 °C±0.5 °C. An 
ambient temperature 3 degrees below the chamber 
temperature caused a noticeable change in the curve 
of viscosity versus rate of shear at high rates of shear. 
The zero shear viscosity at this temperature varies 
about 5 percent per degree. 

Periodically, degradation checks were made on 
samples of 10-percent PIB by taking viscosity versus 
rate of shear data. No changes were observed over a 
2 year period. 



4. Behavior at Small Deformations 

It is clear from eq (2.6) that the shear relaxation 
modulus, G{t) , plays an important role in our curve 
fitting scheme. Since stress relaxation measurements 
cannot be carried out using deformations small enough 
to yield an infinitesimal modulus, this is ordinarily 
obtained by some sort of extrapolation procedure. An 
alternate method, adopted here, is to calculate a curve 
from measured values of the dynamic modulus, C (oo) 
[6]. We employed the conversion scheme given by 
Marvin [7], 



G(t)=G'(co) 



\l(o = t 



H(t) 



(* + T 2 



e-« T W/ln r, 



(4.1) 



using an approximation technique similar to the famil- 
iar derivative schemes [6| for the relaxation spectrum, 
H(t). Even though we must depend on extrapolation 
for values of H{r) at small r, the effect on G(t) is slight. 
Figure 1 shows the dynamic data and figure 2 the 
derived values of G(t). The dotted portion of figure 2 
represents values extrapolated to times smaller than 



4 Certain commercial equipment, instruments, or materials are identified in this paper 
in order to adequately specify the experimental procedure. In no case does such identifica- 
tion imply recommendation or endorsement by the National Bureau of Standards, nor does 
it imply that the material or equipment identified is necessarily the best available for the 
purpose. 
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FIGURE L: Dynamic modulus G'(co) and dynamic viscosity 7)' ((d) 
for a 10-percent PIB (I -100) solution in cetane. 

Solid circles arc the calculated values from eqs (4.2) and (4.3). 




where r\ '(co) is the dynamic viscosity. The solid circles 
in figure 1 represent these values, which agree with 
the measured curve. For values of oo higher than 10 s~ l 
the integrals in (4.2) and (4.3) depend strongly on the 
extrapolated region of G(t). Some of the values of 
G{t) calculated from (4.1) are given in table 1. 

TABLE 1. Shear relaxation modulus for 10 percent 
Vistanex L-100 (in cetane) at 25 °C 



t 


G(t) 


/ 


G{t) 


(s) 


(W 


(s) 


(N/m<) 


1 x 10" 5 


2.11 XlO 3 


3 X10" 1 


3.85 X 10' 


2x 10- 5 


1.98 X10 3 


5 x !()-» 


2.02 X 10' 


5x lO" 5 


1.78 x L03. 


8 X 10-' 


1.05 x 10' 


1 X 10 4 


1.61 x I0 3 


1.0 


7.8 


2x 10-' 


I.45X I0 :5 


2.0 


2.57 


5x 10 -» 


I.I6X LO 3 


3.0 


1.17 


1 X 10 3 


9.5 x 102 


1.0 


4.1 x 10 ' 


2x lO" 3 


7.4 x 102 


8.0 


7.0 x 10 2 


5X 10 - :{ 


5.2 X 10-' 


10.0 


4.0 x LO- 2 


1 x LO- 2 


3.8 x LO 2 


1,1 X 10' 


1.2 x LO" 2 


2x 10" 2 


2.85 x LO 2 


1.8X 10' 


3.9 x lO" 3 


5x 10 - 


1.8 x LO 2 


2. Ox 10' 


2.0 x LO- 3 


1 X 10 ' 


I.08X 10 2 


3.0 X 10' 


9.87 x 10 s 


2X 10 > 


5.92 X 10' 







5. Single Step Stress Relaxation in Shear 

In section 2 the definition of a single step stress 
relaxation history in simple shear was given in which 
the step is instantaneous. In practice, it takes a finite 
time to reach the desired deformation, and it can be 
shown that the measured stress is an upper bound of 
the stress relaxation function. At long times the differ- 
ence between the relaxation function and the bound 
becomes small. With a knowledge of the actual strain 
history, one can estimate this difference and obtain a 
better bound for the relaxation function. 

With the history 



?(r)=0, 







y(r)=fT, 

1 1 



O^r^i (5.1) 

y(r)=y(t) = y, /, ^r^ t 



Figure 2: Shear relaxation modulus. G(t), for a 10-percent PIB where t\ is the time required to reach the maximum 

(L-100) solution in cetane. gtrain ^ eq (2 J) becomes 



The dotted portion shows the extrapolated region. 



o-(/) 



(1/ oj) for the highest measured frequency. As a critical ^' 

check on both the conversion and the extrapolation, 

we calculated C'(o>) and t\'(oj) from the well known _j_ I jy I t\ ~^ + £ >A ( j^t h Z j \ it fc 

relations [6] ) t -t, *\ h )\ h I 



2) 



C'(a>) = cy I tG 
t)'(o>)=J tG 



(I) sin go/ // In / (4.2) Since the quantity given by the integral in (5.2) is nega 

o-(t). 



(t) 



tive, is an upper bound of W(y, t). Since — ^(y, /) 

COS totd In/. (4.3) y . • r • r i • r 

is a monotone nonincreasing function ol /, the ratio of 
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<r(t) 



IW(y, t) goes to unity at high values of t. 



For values of y such that 







] 



a condition satisfied in the stress relaxation data re- 
ported here, it can be shown [8] that cr(t + ti)ly for 
t ^ t\ is a lower bound of W(y, t). 5 

An example of the two bounds and the function 
W{y r , t) is shown in figure 3 for 7 = 0.2 and t\ = 0.03 s. 
The values of W(y, t) are the values that we used to 
correlate our different experiments. We see that even 
for values of £=10 ^i the upper bound is high. The 
average of the two bounds gives agreement to within 
1 percent. In the example in figure 3, a very good 
approximation of W{y , t) can be obtained by averaging 
the time of the two bounds, viz, 



<r(t + iti) 

7 



W(y,t). 



In our experiments the strains were induced through 
a shaft of the rheogoniometer normally used for sinus- 
oidal deformation histories. The motion of the cone 
was monitored with a transducer. The deformations 
were obtained by using a springloaded lever arm to 
drive the shaft through part of a cycle. This arrange- 
ment allowed us to obtain strains up to 7=1.8. By 
using only a part of the sinusoidal deformation, we 
could obtain a motion of the cone very close to a ramp 
function. The time required to reach the maximum 
strain was of the order of 0.01 to 0.05 s. 




LOG ( t,s) 



Figure 3: Bounds ofW(y, t) for y = 0.2 tm<2ti=0.03 s. 



Unfortunately, in our experimental system there is 
another complication due to the motion of the upper 
platen. The true strain at any time t is more nearly 



g(t)- 



K 



where K is a constant depending on the 



stiffness of the torsion bar and the geometry of the 
cone and plate, and g(t) is the nominal strain at t cal- 
culated from the motion of the cone. With a stiffer 
torsion bar the error due to the motion of the upper 
plate would be smaller but small stresses could not be 
measured with enough precision. If the motion of the 
cone is monitored, one can calculate the error due to 
the motion of the plate. In our measurements the error 
in the stress at the early times was 6 percent. 

Figure 4 shows a plot of W(y, t) versus y in which 
both corrections were made. The isochrones at the 
early times are in very good agreement with eq (2.6) 
which corresponds with the lines. To the extent that 
the isochrones are parallel in this type of a plot, i.e., 
superposition occurs by a vertical translation, one can 
justify the representation of W(y, t) as a product of 
the stress relaxation modulus and a function of shear, 
7, i.e., eq (2.6). 



*3 

o 

P , 



t=0.06S_ 
0.08 — 
O.I — 



0.2 



0.3 — 
0.4 — 



0.6 
0.8 
1.0 

1.5 

2.0 




-2 



•This is equivalent to taking the relaxation from the instant of reaching the final shear 
value. 



LOG JJ 
FIGURE 4: Isochrones o/'W(y, t) versus y. 

The lines indicate the calculated behavior from eq (2.6). 



6. Suddenly Applied Constant Rate of Shear 

With our experimental arrangement we found that 
we could not obtain an instantaneous constant rate of 
shear. For that reason we monitored and measured 
the motion of the cone with a transducer during the 
early part of the shearing. The monitoring arrangement 
prevented the use of the temperature control chamber. 
Therefore, the measurements were done at room tern- 
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LOG (t, s) 



FIGURE 5: The so/id lines represent (he stress calculated for a suddenly applied rate 
of shear for the rates of shear indicated. 

The open circles are experimental data, and the solid circles arc calculated values using the true history of the 



perature. Since it took about 20 s to complete a run, 
the temperature variation was no more than ±0.5 °C 
at 25 °C. 

For these runs, due to the limitations of the trans- 
ducer, we could monitor strains only up to y= 0.4. On 
repeating an experiment, different strain histories were 
invariably obtained (over the early part of the time 
scale); this was probably caused by differing positions 
of the gear teeth. So, from each series of repetitive 
runs, we selected the one with the least deviation from 
a constant rate of shear. 

Using eq (2.1) with the history 



y(r)=0, t<0 



y(r) = yr,T^0; y = constant 



(6.1) 



we obtain 8 



<T{t) 



7 1 PWf.fl 



1 + 



d In W(y£, g) 
d In y 



<% (6.2) 



where y= y£ and y is the rate of shear. In figure 5 we 
show the data as obtained. The lines represent the 
prediction of (6.2) and (2.6) which assume a suddenly 
applied steady shear (for the rates of shear indicated), 
and the black solid circles are values calculated from 
(2.1) and (2.6) using the true history of the motion. The 
agreement is very good. The inertia of the upper platen 
and assembly was not corrected for. We estimated for 
y = 22.2 s _1 that the inertial effect caused an error of 



less than 2 percent. Figure 5 shows at high rates of 
shear an overshoot of the stress tr(t) as it approaches 
its limiting value. 

From eq (6.2) the maximum value of the stress should 
occur when 

4tB^--l. ,6.3, 

d In y 



Our empirical form of W(y, t), eq (2.6), satisfies this 
condition at y = 2.44. 

In table 2 we show the strains at which the maximum 
of the stress occurs from other data. 7 Middleman's 
data [9| were estimated from the published graphs. The 
data of the more concentrated polyisobutylene solu- 
tions were obtained in our laboratory. The zero shear 
viscosities at 25 °C were 4460 poise for the 15.1 
percent concentration and 17,760 poise for the 19.3- 
percent concentration. Assuming the principle of 
reduced variables with respect to concentration to be 
valid [10J, we see that the rates of shear reported in 
table 2 for the higher concentrations will be 8 times as 
high when reduced to 10-percent concentration, for 
which the zero shear viscosity is 540 poise. 8 

The position of this overshoot can be used to study 
the error caused by the motion of the upper platen. 
It is easy to show this error by experiments with torsion 
bars of different diameters. In figure 6 we show the 
data obtained from a 15.1 percent solution at a nominal 



7 If eq (2.6) is correct for all concentrations, the i 
the same value of y. 

"The reduced rate of shear y« is given as 



nun of the stress must occur at 



ft j n ^ 

"Throughout the paper we shall use the notation — : to represent the derivative of 

r) In y 
In W with respect to the derivative of the logarithm of the first argument. 



where y\\ is the zero shear viscosity at concentration C\. 
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LOG (t,s) 



FIGURE 6: Stress versus time, using different torsion bars for the some 
nominal rate of shear of 22.2 s~ l applied at t = 0, for a 15.1- 
percent solution of PIB in cetane. 



rate of shear of 22.2 s _1 . The open circles show the 
data obtained with a 0.125-in torsion bar. The solid 
circles are data obtained with a 0.25-in bar and the 
black triangles, with a 0.35-in bar. The histories of the 
motion of the cone for all the runs were almost identi- 
cal, so that a direct comparison can be made of the 
effect of the motion of the upper platen. The maximum 
occurs at t = 0.11 s for the 0.25-in and 0.35-in bars and 
at £ = 0.16 s for the 0.125-in bar. 



TABLE 2. Values of y at which the maximum of the 
shearing stress occurs for the designated rates of shear 





15.1 percent 


19.3 percent 


Middleman's 


y 


Vistanex L-100 


Vistanex L-100 


Data on Polyiso- 




(in cetane) 


(in cetane) 


butylene BASF 




at 25 °C 


at 25 °C 


B-50 at 23.9 °C 


0.0048 






2.45 


.012 






2.4 


.024 






2.3 


1.77 


2..55 


2.55 




5.56 


2.55 


2.55 




11.1 


2.5 


2.5 




22.2 


2.5 


2.55 





7. Steady State Shearing Flow 

In the previous section we examined the stress from 
suddenly applied constant rate of shear. After a certain 
time, as can be observed, the shearing stress reaches a 
limiting value independent of time. The viscosity, 
77(7), which is the ratio of this shearing stress to the 
rate of shear, is given by 



y(y) 



■/■ 

Jo 



PW£, i) l + 



d I n WiM, J) 
d In y 



d£. 



(7.1) 



In figure 7 we show the experimental values of viscosity 
as the open circles. The solid circles are the values 
calculated from (7.1) and (2.6). In the same figure we 
show the dynamic viscosity r/'(co). We observe that 
the two curves cannot be superposed by a shift along 
the horizontal axis. 



8. Stress Relaxation After Steady Shear 

In this history the material is initially sheared at a 
constant rate of shear y until the shearing stress is 
independent of time. The motion stops at time we 
call zero, and we observe the subsequent decay of the 
stress as a function of time. The strain history at time 
r is given by 



y(r)=yo, r^O. 

From (2.1) and the conditions of (8.1) we obtain 9 



(8.1) 



o-(f) 



■\: 



W(y{£-t)*€) 1 + 



d In W 
d In y 



di. 



(8.2) 



In figure 8 and figure 9 we show the experimental 
results with the open circles and the calculated results 
with the solid circles. The experimental points fall 
slightly higher than the calculated values, as the 
motion of the upper platen would lead one to expect. 
At the limit of zero rate of shear, (8.2) reduces to 



o-(Q = f« 

y it 



G(£)d{. 



Differentiating (8.3) with respect to t we get 



(8.3) 



(8.4) 



In principle we can obtain G(t) from this limit, but 
this method is not an accurate one. Since in the previ- 
ous section we mentioned the errors due to the motion 
of the upper platen, it is interesting to look at an 
example where these errors are very pronounced. 
From (8.3) using an integration by parts we obtain 

- " a(t)dt= I G(t)tdt. (8.5) 

y Jo Jo 

From eq (4.2) by taking the limit as co goes to zero, 
we get 

^-f"C(,)i* (8.6) 

ir Jo 



lim 

<o-»0 



9 For brevity, when W appears in an integral without its arguments, it is understood that 
they are the same as for the preceding W. 
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LOG (jf, O, s" 1 ) 
Figure 7: The viscosity as a function of rate of shear. 

The open circles are experimental data, and the solid circles arc calculated values. The dynamic viscosity is also 

shown for comparison. 




LOG(t, s) 



LOGO, s) 



FIGURE 8: Stress relaxation after steady shear, for rates of shear 
y = 0.556 5-M/), y=L77 s-HU), and y= 11.1 sHUI). 

The open circles are experimental data, and the solid circles are calculated. 



FIGURE 9: Stress relaxation after steady shear, for rates of shear 
y = 0.177 s- l (IV),y=l.U s~ l {V) t and y = 5.56 s~HVI). 

The open circles are experimental data, and the solid circles are calculated. 



and from (8.5) and (8.6) we obtain 

(r(t)dt= hm : — 



In figure 10 we show a plot of 



i 



y (i(t)dt 



(8.7) 



versus y using different torsion bars. The error intro- 
duced by the motion of the platen is very dramatic. 

9. Discussion 

The good agreement between the experimental 
data and the calculated values shows the validity of 
the BKZ theory even for materials where the shear 
relaxation modulus falls very quickly. The relatively 
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LOG 
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Figure 10. Variati 



ofy [' 

Jo 



o~(t)dt as a function of y using dif- 



ferent torsion bars, where rr(t) is the measured stress after the 
cessation of steady shear. 



simple form that we used to represent the function 
W{y, t) seems to describe the behavior adequately, 
at least within the range of shearing deformations of 
our experiments. The separation of W(y, t) into a 
product of two functions, one of time and one of strain, 
i.e., W(y, t)=G(t)W r n(y), may be an oversimplifica- 
tion. Although a more complicated form of W(y, t) 
gave us better agreement, we felt that the uncertainty 
in our experimental data did not justify its use. 

Our purpose in this paper is not to establish a unique 
form of W(y, t) , but rather to show how a wide 
variety of measurements from different shear histories 
can be related. Indeed, we know that the behavior 
of plasticized PVC is not consistent with a W(y, t) 
which is a product of a function of strain and a function 
of time. Experiments on normal stresses are more 
critically dependent on the values of W(y, t) at long 
times and large strains. For instance, the first normal 
stress difference, crn (t) — azzit), for a suddenly 



applied constant rate of shear, given by eq (6.1), can 
be calculated from [4] 

(9.1) 

In general, if the function W is such that at high values 
of y we have 



d In W 
3 In y 



<-2 



(9.2) 



then the maximum of cru(t) — & 22 (t) will occur at 
time t m for which 



6 In Wjytm, t m ) 



= -2. 



d In y 

For longer times <xn (t) — (Tiiit) will tend to level to 
a lower value. This behavior of the first normal stress 
difference is roughly similar to that of the shear stress. 
There are clear differences however; a comparison of 
eq (9.1) with eq (6.2) shows that the overshoot for 
the shearing stress, if it occurs, must do so at an earlier 
time than t m . It may well be possible, however, to 
find a form of W{y, t) compatible with the shear 
stress data of this paper which never satisfies the 
condition (9.2) for any value of y. 
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